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Abstract 



The problem of the cosmic ray origin is discussed in connection with their acceleration in supernova remnant 
shocks. The diffusive shock acceleration mechanism is reviewed and its potential to accelerate particles to the 
maximum energy of (presumably) galactic cosmic rays (10 18 eV) is considered. It is argued that to reach such 
' energies, a strong magnetic field at scales larger than the particle gyroradius must be created as a result of the 

acceleration process, itself. One specific mechanism suggested here is based on the generation of Alfven wave 
at the gyroradius scale with a subsequent transfer to longer scales via interaction with strong acoustic turbulence 
in the shock precursor. The acoustic turbulence in turn, may be generated by Drury instability or by parametric 

(N ■ 

\q ■ instability of the Alfven waves. 
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1 Introduction 

O 

Cosmic rays (CRs) were discovered almost a century ago, yet their origin is unknown. There is sufficient evidence 
that at least part of their spectrum (E < 10 18 eV) originates in the Galaxy, while particles of higher energies are 

> ■ 

• <"H . thought to come to us from extragalactic sources. One simple theoretical argument is that particles with such high 

energies (usually referred to as ultra high energy CRs, UHECR) could not be confined to the Galaxy, on account 



of their large gyroradius. Second, their (power law- E~ 2J ) spectrum is harder than that of (presumably) galactic 
particles at E < 10 18 eV (E~ 3 A ), which is consistent with their extra-galactic origin. The last point becomes clear 
if we turn to the other break on the overall CR spectrum, the one at E ~ 10 15 eV, commonly known as the "knee", 
shown in Fig[2 The spectrum above this energy steepens (from E~ 2J to 1 ), so that the premise of their 
extragalactic origin would require an explanation of why the galactic part of the spectrum terminates, while the 
extragalactic part appears exactly at the knee energy. As we shall see, to explain the cosmic ray power law spectrum 
between the "knee" atl0 15 eV and the "ankle" at 10 18 eV in terms of acceleration within the Galaxy is one of the 
most serious challenges of contemporary acceleration mechanisms and one of the main motivations of this study. 

The explanation of the spectrum beyond the ankle (the highest energy event observed so far is 3 • 10 20 eV) poses 
a major challenge to fundamental physics. Given the distance of possible accelerators (at least a few tens of Mpc), 
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particles of such high energy should have experienced significant losses through their interaction with the cosmic 
microwave background radiation (the so called Greisen, Zatsepin, Kusmin or GZK cut-off) while propagating over 
long distances. Note that UHECR are sub-atomic particles with the energy of a well-hit baseball. This energy 
exceeds (by at least three orders of magnitude) that achievable by all existing terrestrial accelerators ( e.g., Large 
Hadron Collider, LHC-10 17 eV). Another fundamental aspect of the problem of CR origin is that they are the 
ultimate receptacle of a significant portion of the gravitational energy in the Universe. Indeed, star formation from 
the gravitational collapse of the primordial gas with subsequent SN (supernova) explosions and their blast waves 
results in CR acceleration. On the observational side, there is an evidence |l]|2][3] (in the form of both synchrotron 
and inverse Compton radiation) that electrons of energies up to 1007W are accelerated in the supernova shock 
waves. 

2 Acceleration mechanism 

The leading CR acceleration mechanism, namely the diffusive shock acceleration (DSA, also known as the 1st 
-order Fermi mechanism) was proposed originally by Fermi in Ref.|4|, and in its modern form by a number of 
authors in the late seventies (5J |6| [7J |8) . The mechanism is basically simple - particles gain energy by bouncing 
between converging upstream and downstream regions of the flow near a shock wave such as that from an SNR 
(supernova remnant) shock. This mechanism requires magnetic fields. First, the field binds particles to the acceler- 
ator (shock wave), in the direction perpendicular to the field. Confinement in the direction along the field lines is, 
in turn, ensured by particles themselves through the generation of Alfven waves by accelerated particles streaming 
ahead of the shock. This occurs via Doppler resonance £2 = k(p/m)p, where £2 and p are the (nonrelativistic) 
gyrofrequency and momentum, k is the wave number, m is the particle mass, and p is the cosine of its pitch angle. 
These waves, in turn, scatter particles in pitch angle (at the rate v ~ £l(mc/p) ($B/Bq) 2 , where c is the speed of 
light) back and force so that they can gain energy by repeatedly crossing the shock. Particle self-confinement along 
the field is thus diffusive and the diffusivity K ~ c 2 /v is inversely proportional to the fluctuation energy 8B 2 , as 
the fluctuating field is responsible for pitch-angle scattering. However, the mean field Bq ultimately determines 
the acceleration rate and the particle maximum energy since it sets the work done by the induced electric field 
on the particles. Indeed, E max ~ (e/c)u s BoR s , where u s and R s are the speed and typical size of the shock wave, 
such as the radius of the SNR shock. The fluctuating part, SB, is typically assumed to be saturated at most at the 
level bB ~ Bq, which provides pitch-angle scattering at the rate of gyrofrequency, and thus limits the particle mean 
free path (m.f.p.) along the field to a distance of the order of gyroradius (the so called Bohm diffusion limit). An 
important thing to keep in mind is that, due to the resonance condition kp = const, confinement of particles of 
higher energies requires that longer waves need be excited. 

The most critical test of this mechanism is the requirement that it accelerates galactic CRs to the energy of 10 15 
eV over the life-time of supernova remnant shocks. Even with the above "optimistic" estimates of the turbulence 
level, the mechanism passes this test at best only marginally. If the turbulence level is lower, then the maximum 
energy should be reduced proportionally. There are indeed a number of phenomena which may either reduce the 
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turbulence level |9|, or which can shift the turbulence spectrum (in wave number) away from resonance with the 
high energy particles and therefore cause their losses 1 10 1. 

Another reason for concern about this mechanism, at least in its standard (Bohm limit) version, is its seeming 
inability to explain acceleration of particles beyond 10 15 eV. As was discussed above, the cosmic ray spectrum has 
only a break at this energy, and continues to about 10 18 eV where the extragalactic component is believed to start 
dominating the spectrum. 

One approach to this problem is to generate a fluctuating component 8Z? significantly exceeding the unperturbed 
field Bq II II . Physically, such generation is possible since the free energy source is the pressure gradient of 
accelerated particles, which in turn may reach a significant fraction of the shock ram energy. Specifically, the 
wave energy density (Sfi/fio) 2 may be related to the partial pressure P c of CRs that resonantly drives these waves 
through the relation 1121 

(8B/Bo) 2 ~M A P c /pu 2 s (1) 

where Ma = u s /Va ^ 1 is the Alfven Mach number and pu 2 is the shock ram pressure. Of course, when 8B/B0 
exceeds unity, particle dynamics, and thus their confinement and acceleration rates, are very difficult to assess if the 
turbulence spectrum is sufficiently broad. The numerical studies by 1 1 1 1 showed that at least in the case of an MHD 
(magnetohydrodynamic) description of the background plasma and rather narrow wave (and particle energy) band, 
the amplitude of the principal mode can reach a few times that of the background field. Moreover, the authors 
of Ref . 1 1 3 1 argue that in the case of efficient acceleration, field amplification may be even stronger, reaching a 
mG (10~ 3 Gauss) level from the background of a few /jG (10~ 6 Gauss) ISM field, thus providing acceleration of 
protons up to 10 17 eV in SNRs. 

Recently, the authors of Ref. 1 14 1 approached this problem from a different perspective. They considered a 
Kolmogorov turbulent cascade to small scales assuming the waves are generated by efficiently accelerated particles 
on the long-wave part of the spectrum. They obtained a particle maximum energy similar to that of 1 13 1. 

Apart from the excitation of magnetic fluctuations during acceleration process, there is yet another aspect of 
the CR- magnetic field connection discussed in the literature. Zweibel| 15 1 points out that since CRs were already 
present in young galaxies (observed at high redshifts), magnetic fields of appreciable strength must also have been 
there at that time. She emphasizes, however, that the approximate equipartition between the CR and magnetic 
field energy established by the current epoch in our Galaxy is not required for the acceleration mechanism and 
presumably results from the fact that they both have a common energy source, namely the supernovae. Indeed, as 
we discussed already, the magnetic field strength merely determines the maximum energy of accelerated particles, 
given the time available for acceleration and the size of the accelerator. If the latter are sufficient then the total 
energy of accelerated particles can, and in most of the DSA models does, exceed that of the magnetic field. The 
latter, in turn, remains unchanged, apart from the conventional compression at the shock and the MHD fluctuations 
discussed earlier. 

In this paper we discuss the possibility of a different scenario, in which the magnetic field may absorb a sig- 
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nificant part of the shock energy as a result of the acceleration process, which may in fact be strongly enhanced. 
The mechanism of such enhancement is based on the transfer of magnetic energy to longer scales, which we call 
inverse cascade for short, even though specific mechanisms of such transfer may differ from what is usually un- 
derstood as a cascade in MHD turbulence. This transfer is limited only by some outer scale L out such as the shock 
precursor size K(p max )/u s ~ r g (p max )c/u s ^> r g (p max ). This approach is in contrast to the above discussed models 
HI 3 1 1 1 41 . which operate on the generated magnetic fields with the scale lengths of the order of the Larmor radius 
fg(pmax) of the highest energy particles and smaller. The advantage of the inverse cascade for the acceleration is 
that the turbulent field at the outer scale bB(L m „) = B rms (which necessarily must have long autocorrelation time) 
can be obviously regarded as an "ambient field" for accelerated particles of all energies. If B nm ^> Bq, then the 
acceleration can be enhanced by a factor B rms /Bo- Note that the resonance field $B(r g ) may remain smaller than 
B r ,ns, so that standard arguments about Bohm diffusion apply, and it is less likely that this field will be rapidly 
dissipated by nonlinear processes, such as induced scattering on thermal protons, not included in the enhanced 
acceleration model 11 1 111 31 . 

As it should be clear from the above, an adequate description of the acceleration mechanism must include both 
particle and wave dynamics on an equal footing. In fact the situation is even more difficult, since the acceleration 
process turns out to be so efficient that the pressure of accelerated particles markedly modifies the structure of the 
shock (both the overall shock compression and the flow profile). 

3 Accelerated Particles and plasma flow near the shock front 

The transport and acceleration of high energy particles (CRs) near a CR modified shock is usually described by 
the diffusion-convection equation. It is convenient to use a distribution function f(p) normalized to p 2 dp. 



Here x is directed along the shock normal, which for simplicity, is assumed to be the direction of the ambient 
magnetic field. The two quantities that control the acceleration process are the flow profile U(x) and the parti- 
cle diffusivity K(x,p). The first one is coupled to the particle distribution / through the equations of mass and 
momentum conservation 




dt dx dx dx 3 dx dp 



(2) 



dt ^ dx ^ 



(3) 




(4) 



where 




(5) 
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is the pressure of the CR gas and P„ is the thermal gas pressure. The lower boundary /?,■„,■ in momentum space 
separates CRs from the thermal plasma that is not directly involved in this formalism but rather enters the equations 
through the magnitude of g at p = pi„j, which specifies the injection rate of thermal plasma into the acceleration 
process. The particle momentum p is normalized to mc. The spatial diffusivity K, induced by pitch angle scattering, 
prevents particle streaming away from the shock, thus facilitating acceleration by ensuring the particle completes 
several shock crossings. 

The system (I2I5> indicate a marked departure from the test particle theory. Perhaps the most striking result of 
the nonlinear treatment is the bifurcation of shock structure (in particular shock compression ratio) in the parameter 
space formed by the injection rate, shock Mach number and particle maximum momentum 1 16 1. 

4 Wave dynamics in the CR shock precursor 

The transformation of magnetic energy to longer scales, while bearing certain characteristics of the conventional 
turbulent dynamo problem, is still rather different from it, in its conventional form. First, this process should 
take place in the strongly compressible fluid near the shock. Second, the Alfven wave turbulence is generated by 
accelerated particles via Cerenkov emission, and thus is strongly coupled to them. Third, the shock precursor itself 
is unstable to emission of acoustic waves. The latter phenomenon is known as the Drury instability and will be 
discussed later. Acoustic waves, in turn, interact with particle generated Alfvenic turbulence, stimulating decay 
instability ( i.e., "inverse cascade"). 

The spatial structure of an efficiently accelerating shock, i.e., the shock that transforms a significant part of 
its energy into accelerated particles, is very different from that of the ordinary shock, Fig|2] Most of the shock 
structure consists of a precursor formed by accelerated CRs diffusing ahead of the shock. If the CR diffusivitity 
K-(p) depends linearly on particle momentum p (as in the Bohm diffusion case), then, at least well inside the 
precursor, the velocity profile U(x) is approximately a linear function of x, where x points along the shock normal 
1171 . Ahead of the shock precursor, the flow velocity tends to its upstream value, Ui, while on the downstream 
side it undergoes a conventional plasma shock transition to its downstream value U2 (all velocities are taken in the 
shock frame). This extended CR precursor (of the size Lqr ~ K(p max ) /U\) is the place where we expect turbulence 
is generated by the CR streaming instability and where it cascades to longer wavelengths. 

4.1 Alfvenic turbulence 

The growth rate of the ion-cyclotron instability is positive for the Alfven waves traveling in the CR streaming 
direction i.e., upstream, and it is negative for oppositely propagating waves. The wave kinetic equation for both 
types of waves can be written in the form 



dNf 3co ± dNf Beo* dNj* 



dt dk dx dx dk 



(6) 
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Here Nt 1 denotes the number of quanta propagating in the upstream and downstream directions, respectively. Also, 
CO* are their frequencies, co ± = kU ± Wa, where Va is the Alfven velocity. The linear growth rates y* are nonzero 
only in the resonant part of the spectrum, kr g (p max ) > 1, In the most general case, the last term on the r.h.s. 
represents nonlinear interaction of different types of quanta and NZ and, if compressibility effects are present, 
also interactions between the same type. As seen from this equation, the coefficients in the wave transport part of 
this equation (l.h.s.) depend on the parameters of the medium through U and Va, which in turn, may be subjected 
to perturbations. This usually results in parametric phenomena |18|. We will concentrate on the acoustic type 
perturbations (which may be induced by Drury instability), so that we can write for the density p and velocity U 

p = p + p; U = U + U 
The variation of the Alfven velocity Va = Va — Vao 

V A ^- l -V A — . 
2 po 

For simplicity, we assume that the plasma p < 1 (which is not universally true in the shock environment) and neglect 
the variation of U compared to that of Va in eq.0. The above perturbations of Va in turn induce perturbations of 
Nj^, so we can write 

Our goal is to obtain an evolution equation for the averaged number of plasmons (A^). Averaging eq.Q we have 

Here the index x stands for the x-derivatives. To calculate the correlator ^ j^-W^i 11 tne l ast equation, we expand 
the r.h.s. of eq.(|6j retaining only the main linear part in N 

fN+ + C± {iV+,JV,- } ~ -Aai£f£ (8) 

The time scale separation between the l.h.s. and r.h.s. of eq.© suggests that to lowest order, the linear growth 
y + rate is approximately balanced by the local nonlinear term C + . Likewise, the linear damping of the backward 
waves y~may be balanced by their nonlinear growth and conversion of the forward waves C . Generally, the Aco^ 
in eq.<0§ is a 2 x 2 matrix operator. If the wave collision term is quadratic in N, then Aco^ ~ yr. This balance can 
be established only for the resonant waves (j^ ^ 0), whereas our primary focus will be on the extended longwave 
interval k < l/r s (p max ) for which y w 0. In this domain, cascading from the generation region k > \/r g (p max ) 
takes place and the refraction (last) term on the l.h.s. of eq.([7) plays a dominant role along, with the nonlinear term 
on the r.h.s. 

To calculate the refraction term we write eq.©, linearized with respect to jV^, as: 
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where 



L± = li^ u±v ^T x - ku 4k +A ^ 



Solving eq.0 for^V^, from eq.0 we thus have the following equation for (N%) 

|w> + 4w>-»4w>44«>^w> + (c W) > do. 

Here we introduced a diffusion operator for the Alfven waves in k space due to random refraction by the acoustic 
perturbations p (via the density dependence of Va), i.e., 

Dk = \k 2 V A 2 ( —L~ 1 — \ (11) 



4 \ Po Po , 

is an example of the well-known phenomenon of induced diffusion. Transforming to Fourier space, we first 
represent p as 

and note that due to the local Galilean invariance of L, we can calculate its Fourier representation in the reference 
frame moving with the plasma at the speed U (x) as: 

Lf q = ±iqV A + Aco± - kU x ^- (12) 



Then, eq. fll 1> can be re-written as: 

Dk = \k 2 vlY j q 



XL}£ (13) 



The last (wave refraction) term on the r.h.s. of ea.( !12> can be estimated as U 2 /K(p max ), which is the inverse 

Aoxt 



acceleration time and can be neglected as compared to the frequencies qV A and Aco. Hence, for we have: 
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For further convenience, we introduce here the number of phonons 

K = — 
q % 

where W q is the energy density of acoustic waves (with coj = qC s as their frequency) 

Wq s Po 
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For Dk in ea.( !10> we thus finally have 




A< 




2C?Po 



<7 2 V A 2 +Ac^ 



Note that Dk represents the rate at which the wave vector of the Alfven wave random walks due to stochastic re- 
fraction. Of course, such a random walk necessarily must generate larger scales (smaller k), thus in turn facilitating 
the confinement (to the shock) of higher energy particles. Thus, confinement of higher energy particles is a natural 
consequence of Alfven wave refraction in acoustic wave generated density perturbations. 

4.2 Acoustic turbulence 

Unlike the Alfvenic turbulence that originates in the shock precursor from accelerated particles, there are two 
separate sources of acoustic perturbations. One is related to parametric 1 18 1 processes undergone by the Alfven 
waves in the usual form of a decay of an Alfven wave into another Alfven wave and an acoustic wave. The other 
source is the pressure gradient of CRs, which directly drives instability. The latter leads to emission of sound 
waves due to the Drury instability. By analogy with eq.Q we can write the following wave kinetic equation for 
the acoustic waves: 



instability of Alfven waves. Note, however, that the combination of Drury instability and decay instability can 
lead to generation of mesoscale fields at a faster than - exponential rate, by coupling together the Drury and decay 
instability processes. 

4.2.1 Decay Instability 

The mechanism of this instability is the growth of the density (acoustic) perturbations due to the action of the 
ponderomotive force from the Alfven waves. This force can be regarded as a radiative pressure term appearing in 
the hydrodynamic equation of motion for the sound waves (written below in the comoving plasma frame) 




Here jo is the Drury instability growth rate and "ft is that of the decay instability. We first consider the decay 




Eliminating velocity by making use of continuity equation 




we obtain 




(14) 



The Alfven wave pressure can be expressed through their energy 
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k 

Using the relation (|9jl between the density perturbations and the Alfven waves and separating forward and backward 
propagating sound waves p ± , we can obtain from eq. (!14> the following dispersion relation for the acoustic branch 

or on writing co = ±qC s + ij^, we have the following growth rate of acoustic perturbations 

Note that the instability requires an inverted population of Alfven quanta. As they are generated by high energy 
resonant particles in a finite domain of k space, such an inversion clearly can occur. 

4.2.2 Drury Instability 

This instability also leads to efficient generation of acoustic waves and it is driven by the pressure gradient of the 
CRs in the shock precursor. The growth rate has been calculated in Ref.| 19 1 (see also II20II21II221 '). and can be 
written in the form: 

^ = _^ ± W, fc\ (15) 

pK C s p V 3lnp / 

Here Pq and Pq x are the CR pressure and its derivative, respectively, and jc is their adiabatic index. For an 
efficiently accelerating shock jc ~ 4/3. Note that we have omitted a term —U x which is related to a simple 
compression of wave number density by the flow and should be generally incorporated into the r.h.s. of eq. (!14t 
(see 1 19|). This term is smaller by a factor C s /U than the second (destabilizing) term. The first term is damping 
caused by CR diffusion, calculated earlier by Ptuskin |23 1. 

5 Mechanisms of transfer of magnetic energy to larger scales 

As it follows from the above considerations, there are a variety of nonlinear processes that can lead to the transfer 
of magnetic energy (generated by accelerated particles in form of the resonant Alfven waves) to longer scales. 
First, as it can be seen from ea.( llO> (last term on the l.h.s.) scattering of the Alfven waves in k -space due to 
acoustic perturbations transfers magnetic fluctuations away from the resonant excitation region to smaller (but also 
to larger) k. Second, the nonlinear interaction of Alfven waves and magnetosonic waves represented by the wave 
collision term on the r.h.s. can be responsible for such process. It is also well known that in the presence of nonzero 
magnetic helicity there is a strong inverse cascade of magnetic energy |24 25 1. Finally, even in the frame work of 
the weak turbulence, induced scattering of Alfven waves on thermal protons leads to a systematic decrease in the 
energy of quanta which, given the dispersion law, means again transformation to longer waves. 
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Returning to the wave refraction process on acoustic perturbations generated by the Dmry instability, it is 
important to emphasize the following. As is seen from the instability growth rate, ea.(!15>. it is proportional to 
the gradient of P c . As the latter should be increased as a result of instability, through a better confinement and 
faster acceleration, this will reinforce the instability, possibly triggering "explosive" growth. This can significantly 
contribute to mechanisms of regulation of P c discussed earlier in 1 16 1. 

6 Conclusions and Discussion 

We have considered a number of possible mechanisms for generation of large scale magnetic field in front of 
strong astrophysical shocks. All these mechanisms are immediate results of the particle acceleration process. Such 
generation is necessary to further accelerate particles well beyond the "knee" energy at 10 15 eV, as is suggested 
by observations of the CR background spectrum and the wide consensus on their SNR shock origin. The fact 
that accelerated CRs constitute an ample reservoir for the turbulence which is required to further accelerate them 
provides a logical basis for our approach. Indeed, the scenario proposed here may be viewed as a self -regulating 
enhanced acceleration process, which ultimately forces the energetic particle pressure gradient to its marginal point 
for Drury instability. 

A new description of the instability of Alfven wave spectrum to acoustic modulations is given. Along with the 
Drury instability this is shown to provide an efficient mechanism for transformation of magnetic energy to longer 
scales. It should be also emphasized that the theoretical analysis of CR acceleration is a challenging problem 
in plasma wave physics, particle kinetics and shock hydrodynamics. It should, and indeed must, include a self- 
consistent description of particle transport, wave dynamics, and shock structure. 
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Figure 1 : Cosmic Ray spectrum 
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Figure 2: Schematic representation of nonlinearly accelerating shock. Flow profile with a gradual deceleration 
upstream is also shown at the top. 
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